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Abstract: In this paper, we study a holographic dual of a confined fermi liquid state by 
putting a charged fluid of fermions in the AdS soliton geometry. This can be regarded 
as a confined analogue of electron stars. Depending on the parameters such as the mass 
and charge of the bulk fermion field, we found three different phase structures when we 
change the values of total charge density at zero temperature. In one of the three cases, 
our confined solution (called soliton star) is always stable and this solution approaches 
to the electron star away from the tip. In both the second and third case, we find a 
confinement /deconfinement phase transition. Moreover, in the third one, there is a strong 
indication that the soliton star decays into an inhomogeneous solution. We also analyze 
the probe fermion equations (in the WKB approximation) in the background of this soliton 
star geometry to confirm the presence of many fermi-surfaces in the system. 
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1. Introduction 



The AdS/CFT has provided us a powerful method to analyze strongly coupled con- 
densed matter systems. There have been remarkable progresses on holographic realiza- 
tions of Fermi surfaces. For example, the fermions in the backgrounds of charged AdS 
black holes |Q, ^, ^, ^ are found to be dual to non- Fermi liquids, while those in electron 
star solutions [^, ^, |9|, 11, ^ (see also |13, 14 1 for a similar calculation in global AdS) 
are dual to analogues of Landau- Fermi liquids. For closely related examples refer also to 
e.g. [|l5|, |l^, 0, |l^, 1£, In these models, the gravity duals are expected to be dual 
to deconfined phases of the Fermi surfaces [^, 22, 23]. The properties of fermi surfaces 



have been usually detected only for their singlet sectors with respect to the SU{N) gauge 
group. Recently, however, new identifications of holographic deconfined non- Fermi liquids 
have been given in [24, 25, where we can confirm the properties of fermi surfaces even 
for non-singlet sectors by looking at the holographic entanglement entropy |27, |2^, |30| 
and specific heat. We gave a simple gauge theory model which is expected to have 0(A^^) 
fermi surfaces in the appendix^ of this paper. 

On the other hand, in realistic condensed matter systems, Fermi surfaces are confined 
in that there are no gapless (emergent) gauge bosons. Therefore it is very interesting to 
construct holographic duals of confined fermi surfaces as initiated in the works by Sachdev 
1 23, A hard wall model with probe fermions has been studied in |31]. The purpose 
of this paper is to present a fully back-reacted example of confined fermi surfaces by 
employing the five dimensional AdS soliton geometry [3^ . We assume a semi-classical 
description of fermi liquids (Thomas-Fermi approximation) in the AdS soliton background. 
We follow the analysis done in the studies of electron stars 0, 11 1, where charged 



fermionic fluids have been considered in a perturbed Lifshitz geometry [34|. Therefore we 
will call our solutions soliton stars in this paper. Also, our model can be regarded as a 
fermionic analogue of the confined holographic superconductors 36]. Notice that as in 
the electron star models, we can calculate the (free) energy from our classical gravity dual, 
while to obtain the specific heat and entanglement entropy we need to perform one loop 
calculations in our gravity dual.^ 

Below we will find that three different phase structures are possible at zero temper- 
ature, depending on the values of two parameters m and /3 of the bulk fermion field, m 
denotes the mass of fermion normalized properly. /3 is a parameter related to the fermion 
charge and number of species of fermions (see appendix |^). Interestingly, in one of the 



""^This is because in our star solutions the (free) energy is the same order as the curvature of gravity as 
required by the Einstein equation. To reahze this, the chemical potential of the fermi liquid gets infinitely 
larger if we take the large A'' limit. Since the specific heat and (entanglement) entropy are reduced by the 
factor of chemical potential compared with the energy, they get much smaller in the large A'^ limit. 
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three cases, we find that the soUton star solutions approach the electron star solutions 
in major parts of spacetime. In both the second and third case, we find a first order 
confinement/deconfinement phase transition. Moreover we found that in the third phase, 
the soliton star solution gets unstable at a certain charge density and should decay into a 
certain new solution which is expected to be inhomogeneous. 

In this paper we also analyze the probe fermion equations in the background of this 
soliton star geometry. In this analysis we work in the WKB approximation, neglecting 
the momentum along the compact spatial direction. We use this analysis to exhibit the 
presence of a large number of sharp fermi-surfaces in the soliton star system. 

This paper is organized as follows. In §Q we present the basic setup of our model. In 
§^, we give perturbative solutions of the soliton star. In §^ we study electron star solutions 
in AdSs. In §|5|, we numerically analyze the soliton star solutions. In §^, we summarize 
the behavior of soliton star solutions and study their phase structure. In §|^, we mention 
the analysis of soliton star without the U{1) charge. In §^ we consider a probe fermion 
in the soliton-star background in the WKB approximation and provide evidence of the 
existence of fermi surface similar to an ordinary landau fermi liquid. In we summarize 
our conclusions and discuss future problems. 

2. The Basic setup 

We start with the Einstein-Maxwell system defined by the action 

Sem = ^ld'x^R+f,- ^^^^^M.) • (2-1) 

Then we couple this system with fermions. We treat the fermions as a fermion perfect 
fluids with the proper equation of state. ^ 

In this setup, without losing generality, we can simply set k = e = L = 1 via the 
rescaling of parameters of the fermi liquids as explained in the appendix ^ and we will do 
so throughout in this paper. 

2.1 The Equations of motion 

The Maxwell equation is given by 

V^F/=j(/), (2.2) 

and the Einstein equations are 

G^. - Qg,u = - [f^xFI - ^^x?^^'^^/..) + , (2.3) 

if) if) 

where j{; and T^i/ are the charge current and the stress tensor due to the fermions. 

^This may be regarded as a fermionic analogue of the holographic construction of confined supercon- 
ductors studied in |§, ||, ||, 0, |l|, ^, 
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Here the coordinates are taken to be {t, z, 9, x, y}, where the 9 coordinate is considered 
to be compact as ~ + 2tt. The spacetime boundary is X and resides at z = 0. 
The metric and gauge field ansatz for the sohton star are taken to be 

ds'^ = -fiz)dt^ + g(z)dz^ + k(z)d9'^ + 4t (dx^ + dy'^) , 

(2.4) 

A = h{z)dt. 

We take the charge current due to the fermions to be of the form 

ju = a{z)Ufj, (2.5) 

where o"(z) is the local charge density of the fermion- fluid and is the local fluid velocity. 
The stress tensor due to the fermions is taken to be that of the ideal fluid form 

tI,1^ = ip{z) + P{^)) u^,u, + p{z)g^,, (2.6) 

where p{z) and p{z) are respectively the local pressure and energy density of the fermions. 
We choose to work in a frame where the velocity is n = {—7=^, 0, 0, 0, 0}. 



It can be shown that the for solving the coupled set of equations ( p.2| ) and (2^) it is 
sufficient to solve the following set of equations. 



p{z)r{z) + p{z)f'{z) - 2y^)a{z)h'{z) + 2f{z)p'{z) = (2.7) 

zf'iz) izk'iz) - w^m (Hz) - .m) ^ ^ _ ^^^^^^^^ _ ^^^^^ ^ ^ ^^^^ 

2f{z)k{z) 2zf{z) 2zg{z) ^^'^^^'^ 2zk{z) z^ ~ ^ ^^'^^ 

f{z) g'{z) k'{z) 3 



2z^f{z)g{z) 2zmzf 2z^g{z)k{z) z^g{z) 
nz)k'{z) , g'{z)k'{z) k"iz) , k'jzf 



4z'^f{z)g{z)k{z) 4z'^g{z)^k{z) 2z^g{z)k{z) Az'^g{z)k{zf 

(tw + - + 1) + ^v7M.(--).(^-) - 2."(.) = 0. (2.11) 

Here ( ^.11| ) follows directly from the Maxwell equation ( |2.2D , while the rest are some linear 
combination of the Einstein equations (|2.3D and the Maxwell equation. Note that (^]^) is 
purely an equation involving the thermodynamic functions characterizing the fermions. 

2.2 The thermodynamic functions of the fermionic fluid 

Let the local chemical potential in the tangent frame at a point be given by 

p{z) = (2.12) 
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Assuming the semi-classical approximation of the fermi liquid in the bulk, as has been done 
for electron stars [0, ^, the energy density and the charge density to be given by 

p{z) = f3 ES{E)dE, (2.13) 

J m 

a{z) = 13 S{E)dE, (2.14) 

J m 

where ^ is a constant which is proportional to the number of fermi surfaces, where S{E) 
is the density of states. The pressure is given by the Gibbs-Duhem relation 

p{z) = -p{z) + f,{z)a{z). (2.15) 

Note that these expressions for pressure, energy density and charge density in (2.13), (2.14) 



and ( p. 15 ) automatically satisfy the equation (2.7). 



In five dimensions we take the density of states for the fermions to be (we assume that 
the fermions see the local flat space). 

S{E) = E{E^ -m^), (2.16) 

m being the mass of the fermion. With these expression we now have to solve the equations 
( |2.8| ), ( p.9p , ( p.lOD and ( |2.11 ) for the four unknown functions in the metric and the gauge 



field. 



3. Perturbative solution for the soliton-star 



If we treat the /3 parameter in ( p.l3[ ) and ( 2.14| ) to be small then we can construct our 



desired soliton star solution perturbatively in (3 about the AdS soliton solution. In this 
section we present such a construction upto linear order in (3. 

We will adopt the following strategy. Inside the star where f3 is non-zero we shall solve 
the equations linearized in f3 about the AdS soliton solution. Everywhere inside the star 
the local chemical potential {p.{z)) will be greater than the mass (m) of the fermion. At the 
radius [zj.) of the star the chemical potential equals the mass of the fermion p.{zr) = m; in 
fact this is the condition that is used to determine the radius of the star. Outside the star 
z > Zr, we have to solve the vacuum Einstein-Maxwell system with (3 being set to zero. 
We will then patch up both the solutions at the radius of the star ensuring continuity of 
the field strengths. 

3.1 Solution inside the star z < Zr 

We can solve for the function g{z) once and for all using the equation ( p.8D after plugging 
in the pressure as described in §^[^ and obtain 

15/(z)3/2 {z {k{z) {2zh'{zf - Af'jz)) + zf'{z)k'{z)) + 4f{z) {k{z) - zk'jz))) 
^^^^ z^k{z) {-l5/3m'^f{z)^h{z) + 10f3m^f{z)h{z)^ + 8/(z)V2 (^^5 _ 45) _ 3/3/j(^)5) ' 

(3.1) 
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Note that this relation is exact in (3. Then we are left with three equations (which follows 
from (|2.9| ) , ( |2.1C| ) , and ( 2.11] ) respectively) after substitution of the thermodynamic function 
according to the prescription of § |2.2| . These are three second order equations for the three 
unknown functions h{z), k{z), and f{z), then there are 6 initial value (boundary) conditions 
to be specified. 

We solve these equations with the following 6 boundary conditions at the tip of the 
soliton (which we take to be unity zt = 1): 

1. There are no corrections to zt (which is set to 1). This is ensured by the fact that 
the corrections to the leading order functional form of k{z) vanishes at z = 1. Near 
the tip k{z) = ki{l- z) + 0(1 - zf. 

2. f{z) does not have any logarithmic divergences at the tip zt = 1. 

3. /(I) = 1. 

4. h[z) does not have any logarithmic divergences at the tip zt = 1. 

5. h{l) = ho. 

6. ki = 45(_i) where ^(-i) is given by g{z) = [1 — z)~^ + 0{1 — z)^. This is required 
in order that there is no conical deficit of 9 at the tip. 

We intend to solve these equation perturbatively in j3 about the AdS soliton geometry. 
The the solution inside the star, upto linear order in /3, is given by 



fU^) = ^ + MSi'{^) + 0{P% (3.2) 



' in 

9rn{^) = ( ^2(l^_^4) ) + (^) + 0(/3'), (3.3) 

hn{z) = 1 (1 - + /3 k^{z) + 0(/32), (3.4) 

hn{z) = ho + P h'il\z) + 0{(3^). (3.5) 

We must remember that this spacetime exists for z < 1. The first order functions are 
complicated and they are explicitly given in appendix ^ 

3.2 The radius of the star Zr 

The radius of the star (z^) is obtained by the following condition 

Kzr) = ^'"^^''^ = m. (3.6) 

V Jin{Zr) 

The radius of the star has the following form 

Zr = ^ + (3.7) 
ho 

where zi^^ is the 0{f3) correction to the star radius and is again given explicitly in appendix 
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3.3 The solution outside the star z < Zr 

Outside the star we have to solve the same set of equations but now with (3 set to zero. In 
this region the g{z) function is expressed in terms of the other functions as follows 

z {k{z) [2zh'{zf - Af{z)) + znz)k'{z)) + 4/(z) {k{z) - zk'jz)) 
'^'^ = 24z^ f{z)k{z) • 

The zeroth order solution in /3 is identical to the solution inside the star. However the first 
order solutions are different. We therefore have 

/o«*(^) = ^ + /3/£l(^) + 0(/3'), 

9out{z) = (—-}—-] + 13 5« (z) + 

(1 )J (3.9) 

1/1 



Kut{z) = h^ + ^h^^l^{z) + 0{^^). 

Here g^^ti^) is expressed in terms of the other functions through the relation ( |3.8| ). The 
rest of the functions are given by 

. (3.0) 
,.(■) ,.,1 _ P2 - 2) - (z* + 1) log (z- - 1) + 24 (A-i (z- + 1) + iK2 (z- - l)) 



h^oltiz) = -\hi log (1 - z") + log [z^ + 1) + H2. (3.12) 



Where the integration constants Fi,F2, Ki, K2, Hi and H2 are determined by condition 
that the functions f{z), g{z), and k{z) and their first derivatives are continuous at the 
star radius Zr- This provides us with 6 conditions to solve for the 6 integration constants. 
Again their explicit expressions are given in appendix |^. 

The behavior of the function for a typical values of parameters is shown in fig.|^. We use 
this perturbation theory as a check on our numerical calculations, which is presented later 
in §^. For the values of values of parameter used in the plot fig.|l] we find good agreement 
with the numerical result as shown in fig.^ 

3.4 Comparison with other solution 

In order to decide which is the dominant solution for a given parameter range, it is necessary 
to consider all existing solution for that parameter range. In this paper we primarily focus 
on the zero temperature solutions and since we fix the temperature (to zero) therefore we 
have to work in the canonical ensemble. Also at zero temperature the energy and free 
energy are same. Hence, we should consider the energies of the zero temperature solutions 
in order the decide which is the dominant solution. The one with lower energy is the more 
stable solution. 
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k(z) h(z) 




Figure 1: Typical perturbative soliton-star (/3 = 0.001, m = 1, Hq = 4.1). The vertical red line 
denotes the radius of the star which is obtained at 2 = 0.2484. The outer solution is denoted by 
the red line while the inner solution is denoted by the blue line. 



Since our system has conformal symmetry, we have to make this comparison of en- 
ergy after fixing the scaling symmetry. We fix the scaling symmetry, by considering only 
solutions with unit radius of the 0-circle at infinity. We now describe this procedure in 
details. 

We have a metric and gauge field of the form 

ds'^ = -f{z)dt^ + g{z)dz^ + k{z)de'^ + z-'^{dx^ + dy^), 
A = h{z)dt. 

Our solutions have the following asymptotic form 

2 

f{z) = fooz'"^ - ruooz^ + -Qz'^ + ... 

1 



g[z) = -IT + 



(3.13) 



,2^--- (3.14) 
k{z^ = kcx^z + z . . . 
h{z) = fi-Qz'^ + ... 

Let us make the following scale transformation on the coordinates 

1111 , , 

z — )• -z; t — )■ -t; X — >• -x; y — -y. (3.15) 

c c c c 

Under this coordinate change the new parameters of the solution are 

/ ^oo ^/ Q I ^ fi f 7/ 2r, ^'1^P.\ 

"^oo ~ 4~' — ^) — loo ~ Joo] K — C Koo- [6.10) 

If we want to set k' = 1 we have to choose 



(3.17) 



oo 
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The time-time component of the boundary stress tensor ( Tu, which we take to be the 
energy of our solutions), should scale in the same way as moo- Hence the scale invariant 
quantities that we should compare is Tuk"^ at fixed the charge density 

Q = -^kL (3.18) 

V /oo 

However, practically it is more convenient to look at the dimensionless ratio 

n^j^, (3.19) 

at constant Q. We therefore compare plots of TZ versus Q for various solutions. The time- 
time component of the energy momentum tensor (i.e. the energy density) for the case of 
out soliton-star is given by 

Tu = ^ + (3.20) 

/oo "'OO 

The chemical potential also provides us a useful information and we will employ the scale 
invariant combination 

/x = -^V^. (3.21) 

V /oo 

We can immediately compute this ratio for extremal RN black holes in AdS. The 
extremal black hole has one parameter less than our soliton star solution and therefore for 
it this ratio TZ is fixed to a definite value independent of Q. The extremal black hole is 
given by the following exact solution 

kiz) = 

/i(^)=^(f"^)-gW^2_ 

with the additional condition f'{z = zr) = 0, zh being the horizon radius defined by 

f{z = zh) = 0. 

For this solution, using standard AdS/CFT prescription, we get = 3m^''^^ and using 
this we get 

7^ = 81. (3.23) 

for all values of Q. 

3.5 Validity of the perturbation theory 

In the perturbative solution for the soliton-star ho may be considered as a parameter of 
the solution. Since /iq is the chemical potential at the tip of the soliton, therefore the 
value of ho must be greater than m for the star to exist. The perturbative solution naively 
exits for all values of ho > m. However as we go on increasing the value of ho the first 
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order correction to various physical quantities increases. Hence, at some value of ho the 
perturbation theory fails to be valid. This point of breakdown of perturbation theory (/ig^^ 
say) may be estimated using the following condition 

/3|z«(/io)|«^. (3.24) 

This condition demands that the first order correction to the radius of the star is always 
less than the leading order value. 

4. The Electron-star 

For comparing the phase structure of various solutions existing in this parameter range 
we must also consider so called electron-star solutions. They have fermionic hairs and 
approach the Lifshitz metric in the IR geometry, which have been first constructed in 
in 4-dimensions ^, ^ (see also recent developments of electron stars in |44, 11, 45, ^]). 



In this section we construct similar solution in 5-dimensions. These solutions have more 
symmetry than the soliton star^ and the function k{z) is fixed to be 

k{z) = 1, (4.1) 
and is no more dynamical. This choice of the k{z) function automatically solves equation 



(2.10) after the thermodynamic functions according to § |2.2| has been substituted. Then we 



are left with three equations to solve for the three unknown functions f{z), g{z) and h{z). 
4.1 IR Lifshitz geometry 

Just like the 4-dimensional case, a Lifshitz geometry emerges in the IR of the electron-star 
which is independently an exact solution of the system of equations. This solution is given 

by 

/(-) = i 

g{z) = f (4.2) 



a 



where gi and hi are given by 



, /a — 1 

hi 



9L 



a 



16Va - la^/^m^ + 5a^ [n? - l) (SOm^ - 53) (4.3) 



60(a(m2 - 1) 1)^ 
o? (-5m^ - 20m2 9) (60m^ - h^w? 3) + 36^ 



^For comparison with our soliton-star we will ultimately consider the ^-direction to be compact which 
breaks this symmetry. However this is not reflected in the metric which retains the symmetry between Q,x 
and y coordinates. 
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where the Lifshitz coefficient a is implicitly given in terms of the lagrangian parameters /3 
and m by the relation 



p = ^720Va - la^/^^ ^ ^Sa"^ + 9q^ + 3a^ - 53a + IGVa - 1 



(4.4) 



A plot of the dynamical exponent a versus (3 for various values of m is shown in fig.^. 
we take a to be greater than 1. Also note that if we assume a to be very large then 
takes the following form 

144 / 1 \ 

/3 



Here 



Remembering that /3 and m are lagrangian parameters which determine a, it is clear 
the above expression that a tends to infinity as (3 tends to zero (from below). This 
corresponds to the extremal black hole solution with an AdS2 throat. Also remember 
the equations of motion assume the relation /x(z) > m, that is, /il > m from ( |4.2| ). It 
a bound < m < 1 for the existence of solutions with this form. 



(4.5) 

from 
limit 
that, 
gives 



200 



150 



100 




Figure 2: A plot of the dynamical exponent a versus /3 for various values of m. 



As emphasized previously this is an exact solution of the Einstein solution. However 
to obtain our electron star solution we shall perturb away from this solution keeping it in 



the deep IR. We shall then use this perturbation result to set the boundary condition for 
our numerical solution. 



4.2 Perturbation away from the IR Lifshitz geometry 

We now consider a perturbation away from the IR Lifshitz geometry in ( |4.2| ) in the following 

way 



9(^) = '^[^ + E9^^'''] (4.6) 



I 



Here /il^ is a parameter of the solution, whose magnitude is unphysical due to a scaling 
symmetry. However, it is very important for it to have a negative sign for the solution to 
exist. Using the scaling symmetry we set its value to -1. 

Also clearly for the perturbation theory to be meaningful we have to consider the 
negative root for 7. The first correction to the Lifshitz solution is given by, 

(1) _ 7 (7 - (r„2 _ 3) + Q, + + 1)^2 + 1) 7) 

^ 72 + 7 + a3 (7 - 3m2) + a'^ (-7 + (7 + 3)m^ + 10) + a (-72 + (72 + 7 - 12) + 19) - 36 
hil^ = ^ - 372 + 97 + a'' (7 - 3to2) + ^3 (-g^2 _^ 3) + 3^2 + (7 - 5)m2 + 3) 

+ a (372 - 67 - 3 (72 - 37 - 4) to2 _ 55) + 36^ ^ ^2(a - 1) (-72 - 7 + (3m2 - 7) 

+a2 (7 - (7 + 3)to2 - 10) + a (72 - (72 + ^ _ 12) - 19) + 36) ^ , 

(4.7) 

where 7 has three roots 
7 = 3 + a, 
1 



+ 2a - Q^m^ - 3am^ - 3 ± 9a'' - 32a^ + 940^ - 128a + 9a^m^ 



2 (a — am^ — 1) 

- 30a^m^ + 25a^m'^ - ISa^m^ + 62a^m2 - 126a'^m^ + 82am^ + 57 J 

(4.8) 

Here 3 + a is the trivial root. As pointed out earlier we have to choose the manifestly 
negative value of 7 (for reasonable values of the parameters) and this is given by 
1 

2 (a — am^ — 1) 



7 =— '-^ ( a^ + 2a - a^m^ - 3am^ - 3 - ga'^ - 32a^ + Ma^ - 128a 



1 

+ 9a^m^ - SOa^m'' + 250^™,^ - 18a^m2 + 62a^m^ - 126a^m^ + 82am^ + 57 " ' 

(4.9) 
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Even in this case we were able to go to higher order but avoid reporting them here to 
avoid clumsy expressions. 



4.3 Numerical solution for the Electron star 



Now we use the analytical perturbative solution about the IR Lifshitz geometry in ( [4 .61 ) 



to set the initial conditions for our numerical solution at point in z in the deep IR. After 
numerically carrying the solution forward we find the radius of the electron star when the 
bulk local chemical potential becomes equal to mass of the fermion. As a check of our 
calculations we have verified that the value of the radius of the star is independent of the 
value of the IR point chosen to put the boundary condition. Also we have to remember 
that there exists an upper bound of the fermion mass which follows from (|4.5| ) and is given 
by < m < 1. 

The solution outside the star radius is again a vacuum solution (with more symmetry 
than in case of the soliton star) and it is possible to write this solution down analytically 



Note that this solution is essentially the black hole solution with an additional parameter 
X which is usually set to unity by a scaling symmetry in the usual black hole solution. 
However in the present case, we have already used this symmetry in the inner solution to 
set /il^ = — 1. The four constants in ( [4.10D , x-, fn^'^'^\ n^^^^ and Q^*^*-* are determined by 
matching with inner solution at the radius of star. 

In Fig]^ we plot the numerical functions inside the star radius for a particular value of 
the parameter. This electron star solution just like the black hole solution does not have 
the extra parameter present in the case of the soliton star. Therefore, like the black hole, 
in this case, the energy density to charge density ratio TZ is determined to be a number 
independent of Q. A plot of IZ versus (5 for various choice of m is shown in fig.§ 

5. Numerical solution for the soliton-star 

Now we go back to the soliton star solutions. To understand its behavior for generic values 
of the parameters, we need to resort to a numerical approach. In order to set the initial 
conditions at the tip of the star for our numerical calculation we first find an analytical 
solution very near the tip of the soliton. 

5.1 Solution near the tip of the soliton 

Let us choose our coordinate scales such that the tip of the soliton occurs at z = 1 just 
as in the case of the perturbative solution in §^. Also similar to the perturbative solution 
we use the invariance under scaling of the time coordinate to set the function /(z) to 1 at 
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Figure 3: Typical solution inside the electron-star for m = 1/3, /3 = 70.43, same values of parame- 
ter for the numerical soliton star plot. Note that the blue line is the actual numerical solution while 
the red line is the perturbative solution about the IR Lifshitz geometry. The numerical solution 
matches very well with the perturbative solution in the IR region but begins to break down when 
z becomes 0(1). This is visible in the plot of the chemical potential \i{z). Here we see that the 
perturbative solution is unable to capture the radius of the electron star and therefore to do so we 
have to use the numerical solution. The yellow line in the /i(z )-plot denotes the value of the mass 
TO = 1/3. 



the tip z = 1. We then make the following ansatz for our dynamical functions near the tip 
z = l: 



,2 



/(z) = l + /i(l-z) + 0(l-z 

9(.z) = ^ + go + 0{l-z), 

1^ — z (^5.1j 

2 , /-o/i -,a2 



k{z) = 4g(_i)(l -z) + k2{l - zY + 0(1 
h{z) = ho + hi{l -z) + 0{l-z 



z 



|2 



Here /iq continues to be a parameter of the solution and has the same interpretation of being 
the value of the local bulk chemical potential at the tip of the soliton. Now substituting 
this ansatz back into the equations and solving them order by order of {1 — z), we can 
determine the unknown constants in (^]^) in terms of ho (and the lagrangian parameters 
which in this case is /? and m). The result of this procedure is given in appendix We 
have actually gone one higher order in (1 — z) in order to put the initial conditions for the 
numerical calculation to greater precision. However since the constants appearing in that 
order are very large and complicated therefore we refrain from reporting them here. 

5.2 Numerical solution 



We use the perturbative solution in §5J_ to set the initial conditions for the soliton star 
near the tip of the soliton. We then numerically evolve the solution forward. Just like in 
the perturbative solution in §1^, we determine the radius of the star (z^) from the following 
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Figure 4: A plot of the energy to charge ratio TZ versus /3 for various choice of m for the electron 
star. 



condition 

fi{zr) = 4= = m. (5.2) 

Then outside the radius of the star we have the vacuum equations (with (3 = 0). Unhke the 
electron star case given by ( |4.10 ), we do not know the analytic solution even outside the 
star. Therefore we again solve these equations numerically with the boundary conditions 
being provided by the inner solution at the radius of the star. This solution is then 
continued to the asymptotic boundary from where the energy density moo, charge density 
Q and the ^-circle radius are read off using a fitting function of the form (3.14). 

A plot of the dynamical functions, for example, for a particular value f3 = 70.43 and 
m = 1/3 of the parameters is presented in Fig:^ This value of /3 is chosen so that the 
Lifshitz solution has the dynamical exponent a = 2. As a check on our numerics we 
compare our numerical solution with the perturbative solution plotted in fig.|^ in §^ for 
small value of /3 (/? = 0.001). The comparison plots are shown in fig.^ which shows that 
there is considerable agreement between the results. 



5.3 Bounds on /iq 

As pointed out earlier the parameter ho must be greater than the mass m of the constituent 



f(z) g(z) 




Figure 5: Typical numerical soliton-star (m — 1/3, /3 = 70.43, /iq — 0.7). The blue line represents 
the radius of the star at z = 0.275. 

fermions. This is because /iq is simply the local chemical potential at the tip of the soliton 
and for any fermion condensation to take place it has to be greater than the mass of the 
constituent fermions. Therefore /iq clearly has a lower bound for the solution to exist. 

More importantly, it turns out that /iq also has an upper bound. There are two reasons 
for the presence of upper bound. One of them can be understood from the fact that the 
denominator of (7(_i) in blows up at this bound. In fact we loose control over the 

numerical calculation as we approach this point even before we reach it. From (|D.2| ) this 
denominator is given by 

V{ho) = 9^hl - W/Shlm^ - Ifyph^m^ + 16 (/3m^ - 45) . (5.3) 

The upper bound is given by the real root of the polynomial P(/io). This upper bound for 
/iQ is denoted by /ig below. At /iq = the metric gets singular (or divergent) at the tip 
z = 1. For example it is given by /iq = 1.053 for f3 = 70.43 and m = 1/3. 

Another reason of the upper bound is because when f3 is larger than some value (3c 
(this is /3c — 20 for m = 1/3), the solution develops a singularity somewhere between the 
tip and the AdS boundary^ before ho reaches /ig. This upper bound of /iq is denoted by 



*It would be no problem, if there were to exist a point satisfying (5.2), before the singularity. However 
it is not the case here. 




Figure 6: In this plot we compare the numerical solution obtained in this section with the pertur- 
bative solution in /?. The value of parameters that we use are /3 — 0.001, m = 1, /iq = 4.1. The 
dotted lines (green and black) represent the numerical solution while the continuous lines represent 
the perturbative solution. Numerically the value of the radius is obtained to be z = 0.2483 which 
compares well with the perturbative value of the radius z = 0.2484. 

/ig(< /iq). For example it is given by Hq = 0.774 for /? = 70.43 and m = 1/3. Thus for 
(3 > (3c, the true upper bound is given by /iq, while for /3 < /3c it is given by /iq. Note that 
this phenomena only happens for the mass range < m < 1. In other words, for m > 1, 
we always have /3c = oo. 

The behavior of the soliton star solution changes when f3 crosses /3c. If /3 > /3c, as /iq 
approaches Hq, the charge density Q gets infinitely large and the star radius Zr gets closer 
to the AdS boundary z = 0. This means that we can realize a very large star. On the other 
hand, if /3 < /3c, when /iq approaches h^, the charge density remains finite and approaches 
to a certain value Qmax- In this case, the star radius also stays a certain range and does 
not get closer to the AdS boundary. The behavior of the star radius Zr in these limits as 
a function of /3 is plotted at m = 1/3 in Fig.^. It will be an interesting future problem to 
analyze this transition in more detail. 

5.4 Thermodynamical Stability 

It is also helpful to examine the relation between the chemical potential fi and the charge 
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Figure 7: The star radius Zr as a function of /? at m = 1/3. 



density Q. The thermodynamical stabihty requires 

where E is the total energy of our system. If this derivative gets negative, then inhomo- 
geneous configurations where the charge is concentrated on some particular regions can be 
energetically favored over the homogeneous one. As we will show explicit numerical re- 
sults in section 6.1 later, the inequality (5.4) is satisfied at least for a small charge density. 
Generically, when the charge density is larger than a specific value Qc; the condition (^^) 
is violated. For such a value of Q the soliton star solution becomes thermodynamically 
unstable and moreover, a perturbative instability is also well expected in the spirit of the 



Gubser-Mitra conjecture |47, 48 1 

In the case < m < 1, for /3 > /3c, no violation of (|5.4|) is observed, while /3 < /3c, 
there exists a point Q = Qc where it starts to be violated. On the other hand, in the case 
m > 1, for any /3 we find the instability point Qc- We will explain how these instabilities 



affect the phase structure of our system later in section 6.1 



6. Comparison of Energy for various solutions 

As follows from our previous arguments, for a given generic value of (/3,m), we will have 
three different solutions at zero temperature: (i) a soliton star (SS), (ii) an electron star 
(ES)and (iii) an extremal charged black hole (BH), which are all smooth solutions in 
our Einstein-Maxwell-Fermion system. Notice that the electron star only exists when 
< m < 1, while the others exist for any /3 and m. The physically relevant solution at 
zero temperature is the one with the lowest energy or equally the one with the lowest value 
of the ratio ( |3.19| ) for a fixed charge density Q. Notice that extremal charged black 
holes have always the same value 7^ = 81 ( 3.23| ). An electron star takes a fixed value 
which does not depend on Q for a given (/3,m). On the other hand, a soliton star has a 



-18- 



Figure 8: The plots of the chemical potential p, and ratio TZ as functions of the charge density 
Q. The value of m is indicated at the top of each graph. Each graph consists of six colored curves 
which correspond to (from the up to the down) (3 = 0.1 (red), /3 = 1 (blue), f3 = 5 (green), /? = 15 
(yellow), 13 — 40.7 (red), /3 — 70.4 (blue). We also plotted the values of TZ for the electron stars 
with these values of /3 as well as the that for the extremal black hole. 



non-trivial function TZ = TZ{Q), which also depends on (/3,m). As we discussed in section 
5.4 , the behavior of the chemical potential fl, normalized as in (|3.2l| ), is very important to 



know the thermodynamical stability. We plotted the TZ and /i as a function of Q for the 
several values of {f3,m). 

6.1 Phase Structure 



Based on the behaviors of TZ and fl, we can classify the phase structure at zero temperature 
into three cases (simply called case 1, case 2 and case 3) as explained in Fig.^. We present 
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Figure 9: The three possible phase structures for various values of /3 and m. In the right figures 
we show the schematic graph of the ratio 7?, as a function of the charge density Q for the extremal 
charged black hole (BH), the electron star (ES) and the soliton star (SS). We also indicate the point 
where the thermodynamical instability defined by < starts. 



the schematic diagram in Fig. 10 to show which case corresponds to a given value of (/3, m). 
Before we explain each case, we would like to note that the soliton star is the stable solution 
with the lowest energy if the charge density is enough small. This is obvious from the fact 
that the AdS soliton is the lowest energy solution to the vacuum Einstein equation |33[ |. 

The case 1 is characterized by the condition /3 > /3c, where /3c is the lower bound of /? 
which allowed the limit of infinitely large charge density Q — )• oo as discussed in section |5.3| . 
In this case, as we can see from Fig.^, the soliton star solution always has the lowest energy 
(or equally TV). The condition ( ^.4]) is always satisfied and there is no thermodynamical 
instability. Moreover, we notice an intriguing fact that the value of the soliton star 
approaches that of the electron star in the large charge density limit. Actually, as we show 
in the section the solution star solution itself approaches to that of the electron star 
in this limit. This is also consistent with the fact that the case 1 occurs only in the mass 
range < m < 1. 

On the other hand, if /3 < /3c, the system corresponds either case 2 or case 3. In 
these cases, we always encounter the upper bound of charge density Qc above which (^!J) 
is violated and the system gets thermodynamically unstable. The difference between the 
case 2 and 3 comes from the detail of how the instability starts. In the case 2, the instability 
occurs at a enough large value of charge density and the energy density of solition star gets 
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Figure 10: The Schematic phase diagram for various /3 and m. The green region, the white region 
and the blue region correspond to the Case 1, Case 2 and Case 3. Each point with a small colored 
circle and triangle, which is defined by the previous picture, represents an example we performed 
an explicit numerical calculation. 

larger than that of either the electron star (for < m < 1) or the extremal black hole (for 
m > 1) before the onset of this instability. Therefore, in the case 2, there is a first order 
phase transition from the soliton star into either the electron star or the extremal black 
hole. This can be regarded as a deconfinement /confinement transition. 

In the case 3, the instability takes place when the energy of soliton star is still the 
lowest one. In this case, the solition star should decay into a new solution (NS) in order 
to have a sensible phase structure as is clear from Fig.^. It is possible that for a much 
larger value of charge density the electron star or extremal black hole is finally favored. 
This appearance of the new blanch will be related to a perturbative instability of solition 
star solution for a large enough charge density. As suggested from the thermodynamical 
instability, it is natural to expect that the new solution is inhomogeneous with respect to 
the X, y directions. It will be a very important future problem to work out this conjectured 
perturbative instability and new solution. 

In summary, the ground state for each of the three cases changes in the following way 
as we increase the charge density Q: 

Case 1 : Soliton star 

Case 2: Soliton star—)- Electron star (0 < m < 1) or Extremal BH (m > 1) 

Case 3: Soliton star—)- New solution (— )• Electron Star (0 < m < 1) or Extremal BH (m > 1)). 

The phase transition found in the case 2 can be regarded as a confinement/deconfinement 
phase transition. In the case 1, this transition goes smoothly and eventually the solition 
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star degenerates with the electron star in the large charge density limit as we mentioned. 
In the case 3, we first encounter an instability against inhomogeneous perturbations. 



6.2 Electron star as a limit of soliton star 



In this section we explain that the soliton star approaches electron star solution when both 
(3 and the charge density are large, as we already mentioned in section |6.1| . The fact is 
intuitively clear since the bulk chemical potential fi{z) monotonically decreases from the tip 
to the AdS boundary and bounded by ho. If we pack more and more fermions, fj,{z) tends 
to be constant and in that case corresponding solution becomes the Lifshitz one which 
coincides with the IR region of electron star. This also explains the reason why the ratio 
TZ of solition star approaches that of the electron star for large /3 and large Q as observed 
in the previous section. 



In figure 11 we compare the part of the metric g{z) of soliton star solution with nearly 
maximum value of ho and that of electron star when /3 = 70.43 and m = ^. To compare 
both solutions we use rescaling ( p. 16 ) so that radii of both stars are same. Although in the 
region very close to tip, both solutions look different but if we move toward the boundary 
they eventually coincide. 
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Figure 11: Comparison of the solutions in m = 1/3, /? = 70.43, Hq = 0.7744 case. The red line is 
g[z] of soliton star and the blue line is g[z] of electron star. The tip of soliton star is located at 
z = 149.8 



7. Soliton Star without Charge 



It is worth mentioning that Einstein fermionic fluid system without any U{1) gauge field 
also admit soliton star solutions. In this section we consider these solutions. One can 
regard the solutions as the poincare patch versions of the AdS degenerate star solutions 
considered in [14|. In this system, equation of motions are given by just setting h{z) to 
zero in Einstein equations in (2.8)(2.9)( p.l0 ), and fj,{z) appearing at ( p.l3 ),( 2T^ ) is given 
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by 



(7.1) 



where ep \s arbitrary constant (and regarded as Fermi energy at the tip) so that the stress 
tensor of the fluid t/^™'' satisfy V^^/S'*'^ = 0- Alternatively, a U{1) gauge transformation 
results in h[z) = ep, and in that frame ( |7.lD coincides with ( 2.12| ). 



One can numerically solve the equations as in previous case. In figure 12 we plot the 
energy density of the star Tu/ foo as function of e_F in the case of m = 1, /3 = 1. 
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Figure 12: Plot of Tu/ foo as function of ep in the case oi m — I, (3 — 1 



8. Probe fermions in the soliton star geometry 

In this section we shall consider a probe fermion in the background soliton star geometry 
obtained in §^. In this study we shall closely follow the analysis of ^, ^] . We wish to show 
(in a WKB approximation) that there are many well defined fermi-surface momentum for 
our soliton star. Because of the confining nature of the background geometry the gapless 
excitations in the system are only those associated with these fermi-surfaces. Thus our 
system is a very good holographic realization of the standard fermi-liquid. Further using 



a recent result in |49| we argue that Luttinger theorem trivially holds for our system. 



8.1 Presence of Fermi-surfaces from a probe fermion analysis 

The background soliton star geometry consists of a metric of the form 

ds2 ^ -f{z)df + g{z)dz^ + k{z)de'^ + ^ {dx^ + dy^) , (8.1) 
together with a gauge field 

A = h{z)dt. (8.2) 
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Since we want to deal with fermions therefore we make the following choice of the tetrad 
for the metric under consideration 

el = vWy, ef = v^; 4=\/M^; < = e| = -, (8.3) 

z " z 

where the bar signifies tangent space indices. 

In this background geometry we consider a spinor field with the following action 



S = J (fx,f^ i \^r''VM^ - m ipi^j + ... (8.4) 

where the dots in the above action signify appropriate counterterms so as to make the 
boundary theory well defined. Here the covariant derivative is given by 

Dm = Om + ^u^abAiT^^ - iqAM. (8.5) 

Note that here we have used small letters (like a, b) to denote the tangent space indices, 
where as capital letters (like M) to denote bulk coordinate indices. Here, m and q are 
respectively the mass and charge of the probe fermion. Also, we have used 

pa6^[pa^p6]^ (8.6) 



and uJabM here is the spin connection. For the choice of tetrad in ( |8.3D the non-zero 
components of the spin connection is 

/'(^) , k'{z) ,^ 1 , 1 , 
Wf? = --dt\ ujn:. = :au: Uxz = ax; Uyz = ay 

(8.7) 



The Dirac equation that follows from (?^) is simply given by 

(^r''efpM-m)v = 0. (8.8) 

Here are related to the tetrad defined in ( |8.3| ) through the relations 

— ^m9 Vab (o-y) 

and are given by 

v/(^) Vaiz) vHz) 

Now let us make the following ansatz for our spinor field 

^ _ ^ ^i(~uit+pe+kix+k2y)-,j, 21) 

if{z)k{z))-^ 

Here p denotes the momentum of the fermion along the 6 circle. Note that due to the 
anti-periodic boundary condition of the fermions p is an half integer times vr. With this 
redefinition the Dirac equation (|8.8D reduces to 



^ T^dr -m + izK^V ]^ = (8.12) 
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where // runs over only boundary coordinates. In ( 8.12| ) is defined to be the foUowing 
quantity 

K,, = {-—j==l^ + qh{z)^,pMM (8.13) 

Now exploring the rotational symmetry in the xy-plane we set A;2 = and ki = k. 

We would also work in the large charge large mass approximation. Also, for capturing 
a nontrivial dynamics we scale the frequency a; and momentum k so that they are same 
order as the mass and the charge. For convenience we choose a large number 7 and keep 
it explicitly in our equations by considering the following scaling 

m ^ J m, ^7— 5-7, CO ^ 'J CO, k ^ j k. (8-14) 

Note that for simplicity we shall not scale the 6 circle momentum. This would amount to 
neglecting the term proportional to p ( ^.12 ). This assumption physically means that we 
are ignoring any dependence of our spinor field on the 6 coordinate. This assumption is 
also justified from the fact that from the boundary point of view we are interested in the 
2 + 1 dimensional physics that is obtained after the Scherk-Schwarz compactification (on 
the 6-circle) has been performed. 

Once the term proportional to p is dropped assuming large 7, there are only 3 Gamma 
matrices that are involved in ( ^.121) . Following ||ll| we make the following choice of the 
Gamma matrices 

/ n \ _ / n \ / n \ 

(8.15) 






With this choice of the Gamma matrices if we take 

^ = 




then from ( 8.12| ) it follows that the equations for ^'i and decouple and without loss of 
generality we can focus on one of the equations given by 

^ crsdz+^i -m-izka2+ J— (cj + )cri ) Wi = 0. (8.16) 



Note that since equation ( ^.16| ) is a equation for a two component spinor "fi and it is 
identical to equation (9) in |ll|. Therefore although our analysis is very similar to that in 



1 11] but our boundary conditions will be the different which will turn out to be extremely 
crucial. Now since ^'i is a two component spinor let us consider 

^1 = 




Then the two equations that follows from ( p. 16 ) are 



^ --d^^i + 71 - m^i - zk^2 + -^={uj + At)^2 ) = 0, 



'V 77( 



z] 



^.17) 



-25- 



Eliminating <I>2 from the above equations and keeping only leading order terms in 7 (which 
is justified if we are not concerned about the precise energies of the low lying states, see 
1 11] for details) we have 

- 7V(z)$i = 0. (8.18) 
where the Schrodinger potential is given by 

V{z) = 9{z) (m^ + z'k' - ^^^^^j^) ■ (8.19) 

We plot this potential in hg.^ for the solution shown in fig.^ for u = 0. We now proceed 
to perform a WKB study of the equations ( p. 18 ) and ( ^.19[ ). It is clear from the plot in 
fig.|l3| that if A; ^ 0, then the potential is always positive (see the green lines in fig.|l^ for 
example). This regime of momentum and frequency range is uninteresting from the point 
of view of capturing the gapless excitations about the fermi surface. Prom fig. 13 we also 
see that for generically (even away from the a; = considered for these plots) the potential 
has a single turning point which we call z*, where the potential V{z) linearly vanishes. 
Note that for a; = and k = (the blue lines in fig.^3|), is simply the radius of our 
soliton star. There are a small number of cases for which there are two turning points (for 
instance the yellow line for m = 0.33, (3 = 70.43, /iq = 0.7). However, the final result of the 
WKB analysis (existence of fermi surface poles) is same as the case with the single turning 
point. Therefore, we present here the details of the WKB analysis for the case with a single 
turning point. 

For the generic case of a single turning point (at z^) in fig. 13, we see that V{z) < 
when 2; G [z^, 1] while near the boundary of AdS for z G [0,^*], we have V{z) > 0. For 
convenience let us define the following quantities 



X{z)=-f dz,/V{z)- Y{z)=j dz^Z-Viz); X = X{0); Y = Y{1). (8.20) 



Then in terms of X and Y we can write the following most general solution of ( 8.18 ) 

$i(z) = A e^^(") + B e-'^^'\ for z G [z^, 1] 
$i(z) = C e^(^) + V e-^^'\ for z G [z^, 1]. 



.21) 



In order to be consistent with the anti-periodic condition of the fermion about the ^-circle 
we would like to the impose vanishing of the fermionic field at the IR cut off z = l.Thus 
demanding ^i{z)\z=i = 0, yields 

A + B 



A-B 



-itanY (8.22) 



Further we impose ro.cLtch.ing conditions cit the turning point -s^star* 

Using the same proce- 
dure as explicated in appendix A of 0] we find 

C = (^ + i3) cos ^ - i{A - B) sin ^ 

If vr 7r\ (8-23) 

V = - {{A + B)s\ii- + i{A- B) cos - 
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Figure 13: Plot of potential V{z) as a function of z for the soliton star in the a; = limit. We 
plot it for various values of of parameters as indicated in the plots. The blue line represents the 
w = 0, fc = plot in all the cases; therefore the turning point of this blue plot represents the radius 
of the star. We increase the momentum k along the red, yellow and green line (in that order). 

Now we know that the poles of the boundary Greens functions occurs where the function 
vanishes at the boundary. This conditions implies that the poles of boundary Greens 
functions occurs when 

I = ^ (8.24) 



Now using (|]2|) and (|^ in (|]2|) we have the following condition 



vr 



tan ( y + ^ ) + 2e~2^ = 0. (8.25) 



In the large 7 limit this condition implies 

y + ^ = nvr - 2e"^^ ^mr. (8.26) 

Note in that the frequencies solved from this equation do not have any imaginary part 
(which is unlike case II in ||ll[). This is expected since we do not have any horizon in our 
soliton star geometry and therefore there is no dissipation. Also note that in the large 7 
approximation that we are working, the exponential terms are very small and therefore 
may be legitimately neglected. 

We know that in a fermi liquid there are low energy excitations near the fermi-surface 
momentum kp, which implies the existence of a pole in the Green's function at a; = and 
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k = kp. Thus in our case we see there are many fermi surfaces corresponding to fermi 
momentum k^p^ described by the equations 

y(cj = 0,k = k^) = n-K. 

It is particularly interesting to note that the existence of large number of well defined 
fermi-surface is similar to that observed in the gauge theory with adjoint fermions in the 
free limit discussed in appendix 

8.2 Luttinger theorem for the soliton star 

According to the Luttinger theorem ^ the total charge density of the system is 

proportional to the volume enclosed by the fermi-surface times the charge of the fermionic 
species under consideration. Our soliton star lacks the presence of any horizon, and all the 
charge in the system is carried by the fermions themselves. Therefore, Luttinger theorem 



is expected to hold in our system as in the electron star case [11|. 

On the other hand, in the gravity dual, the following relation has been found ^ 
(in d + 1 dimensional bulk space-time) 



i^VF = p-A, (8.27) 



(2vr) 

where q is the charge of the fermion, p is the total charge density of the system, and Vp is 
the total volume enclosed by all the fermi-surfaces. Here A is the anomalous term which, 
if non-zero for a system, implies that the Luttinger theorem is invalid in that system. In 
terms of the bulk fields it is given by 

1 



A2 

where A is the gauge coupling and zin is the end point of the IR geometry. For example in 
the case of a black hole, zjji will be the horizon radius, while for the electron star of §^ it is 
zjR = DO. This relation ( ^.27 ) has been proved under very general assumptions (which are 



valid for our soliton star system) in [^] (see also |12] for an analysis where both confined 
and deconfined fermions coexist). 

For our case zjji = 1 and corresponds to the IR tip of the soliton. From the behavior 
of the solution star near the IR tip (discussed in §|^, it follows that the zt-component of 
the gauge field strength F^* (or the z-component of the bulk electric field) vanishes since 
the time component of the gauge field goes to a constant near the tip of the soliton star. 
Thus with ^ = in this soliton star system, Luttinger theorem clearly holds. 



9. Conclusions and Discussions 

In this paper we studied a holographic model of confined fermi liquid by putting fermions 
in the AdS soliton geometry. We numerically solved this back-reacted problem and studied 
its phase structure when we change the charge density at zero temperature. We also 
constructed perturbative solutions analytically and confirmed that they match with the 
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numerical ones. As we summarized in Fig.^ and Fig.^ in section 6A there are three 
different cases depending on the values of the parameters m and /3. 

The case 1 corresponds to the large /3 and small m. This corresponds to a system with 
a lot of fermi surfaces and the back-reaction of fermi liquids in the gravity dual can be 
kept small compared with other two cases. Therefore we have stable soliton star solutions 
even for large values of the charge density. We made an interesting observation that these 
solutions approach those of electron stars away from the tip. The case 2 is realized when 
both /3 and m are small. In this case, due to the back-reaction of fermi liquids in the bulk 
gravity, the soliton star gets unstable when the charge density is very large. It decays into 
the extremal charged black holes or the electron stars via a first order phase transition. 
For large values of m, the soliton star gets unstable at a relatively small charge density and 
the numerical calculations of the energy and charge density predicts that the soliton star 
will decay into a new solutions instead of the extremal charge black hole or the electron 
stars. This is the case 3. The new solutions are expected to be inhomogeneous. It will be 
a very interesting future problem to construct such new inhomogeneous solutions and to 
understand the nature of this phase transition. 

We also performed a probe fermion analysis and show that there are many fermi 
surfaces, which are expected to behave like Landau's fermi liquids. The understanding of 
all of these properties in the gravity dual from the dual CFT viewpoint is an intriguing 
future problem. Also it will be another interesting future direction to see how these phase 
structures change at finite temperature. 
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A. 0{N^) Fermi Surfaces from SU{N) Yang-Mills 

Let us consider an adjoint fermion coupled to a SU{N) gauge field in 4 dimensions with 
a global U{1) chiral symmetry. In order to find an evidence that this system have 0{N'^) 
Fermi surfaces, we wish to evaluate the partition function of this system and hence compute 
the specific heat in the free limit. Possibilities that gauge non-singlet sectors of fermions 
make fermi surfaces in addition to the gauge singlet one have been already discussed in 
1 22] in different setups. 
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We shall closely follow the basic set up in 53]. The grand canonical partition 
function is given by 



Z{/3,^l) = Tr^exp{-^3{H-^iQ))j, (A.l) 

Q being the chiral charge operator which is a conserved U (1) global charge for this system. 
Here /3 = l/T and /u are respectively the inverse temperature and the chemical potential 
for the U{1) global symmetry. Let us consider this theory on x with R being the 
radius and /3 being the radius. Ultimately we shall be interested in the large R limit 
in which the sphere tends to flat i?^. Here we shall set R = 1, so this flat space limit will 
correspond to /3 <^ 1, that is when the radius of the time circle is very small compared to 
the radius. 

The partition function may also be written in terms of a functional path integral 

Z = f VA^V^e-^^^^^^\ (A.2) 



We have to impose periodic boundary condition for the gauge field and anti-periodic bound- 
ary conditions for the fermion along the time circle. In the free limit of the theory all the 
heavy modes of the fields may be integrated out (by evaluating one loop vacuum diagrams) . 
However, the zero mode of the time component of the gauge field Aq which is strongly cou- 
pled even in the free limit has to be treated exactly. Thus, the process of integrating out 
the heavy modes yields a quantum effective action for this zero mode which may be written 
in terms a matrix integral [52, 53] which has the following form 

Z = i dUe'^-ff^^\ (A.3) 



with U being a SU{N) matrix. The effective action Seff{U) may also be evaluated in the 
free limit by explicitly counting the gauge invariant states in the free theory on x and 
then projecting onto the singlet sector so as to ensure Gauss law constraint. This effective 
action is evaluated to be ]^, ^] 

SeffiU) = - f; - {z,{xn + (-l)"+izy(x", fi)) (Tr{U-)Tr{U^ ") - iV (A.4) 

where x = and z^, Zf are respectively the single particle partition functions for the 
gauge field and the fermionic field, which are given by 

- 2a.3 

- (l_x)3 ' 

3 . . (A.5) 



/ X 2x2 / „ 



Note that for concreteness we have considered a chiral fermion in writing down this single 
particle partition function Zf{x,fj,). However the final result of this appendix also holds for 
a Dirac fermion. 
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Now we intend to perform the integral in (|A.3D by rewriting it in terms of the eigen- 
values (e*^) of the matrix U . The integration measure dU is given by the Haar measure 
of SU{N). Introducing a suitably normalized eigenvalue distribution function p{9), the 
effective action in (A^) may be written as [52, |53| 



n=l 



where pn = JZ.^d9p{6) cos n9, are the Fourier coefficients of the eigenvalue distribution. 
Here Vn is given by 

Vn = -(l- - {-ir+'zf{x^,p)) . (A.7) 



n 

Note that until now the effective matrix integral that we have found is exact in N. 
But now we wish to take a large N limit, and evaluate the integral ( |A.3| ) in the saddle 



point approximation. In this limit there are essentially two saddle points. In ( A.6 ), when 
the Vn are all positive they act as a repulsive potential for the eigenvalues, while when 
they are negative they act as attractive potential. As a result, for positive Vn, a saddle 
point is obtained where the eigenvalues are uniformly distributed, that is pn>i = 0. This 
corresponds to the confined phase. On the other hand, when Vn is no longer positive the 
action is minimized by a tightly clustered configuration of eigenvalues and in the large 
limit this is given by the saddle point = for all n. This corresponds to the deconfined 
phase. 

As we have emphasized for the confined phase to exist we must have Vn > for all n. 
This implies that Vn = will provide us with the line of phase transition in the /xT— plane. 
However, since and Zf are monotonically increasing functions with < x < 1, the most 
stringent condition is obtained when n = 1. Therefore the equation of the curve separating 
the confined and deconfined phases in the ^T— plane is given by the equation 

Zv{x) + Zf{x, p) = 1. (A. 8) 



This curve has been plotted in fig. 14 



In this discussion we would particularly like to focus on the deconfined phase in the 
limit ^ ^ 1, /3 ^ 1 so that /3/u S> 1. In this limit, since the radius of the thermal circle 
becomes small compared to that of the radius (which has been set to 1), we we can 
approximate by flat space. The other condition [Sp ^ 1 ensures that we are in the 
deconfined phase. Since we are in the deconfined phase we have p„ = 1. The effective 



action (A.6) in this phase is then given by, 

Vnipp » 1) = E - - ^-(^") - {-^r^'zf{x\p)Y 

n=l n=l ^ ^ 



s, 



(A.9) 



deconfined 



In the sum in ([a!9|) there are two contributions. One contribution comes from the gauge 
bosons through Zy and the other contribution comes from the fermions through zj. The 
bosonic contribution is independent of p which a consequence of the fact that the gauge 
bosons are uncharged under the global U{1) symmetry. As a consequence of the this fact, in 
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Figure 14: Confinement deconfinement phase transition curve in the /iT-plane 

the limit (/3/x ^ 1) under consideration, the leading contribution comes from the fermions. 
The first term in ( [A.9D merely contributes to a constant additive shift to the effective action 
and may be absorbed in to a normalization of the partition function. 

In the large N limit, the grand potential in this deconfined phase saddle point is given 

by 

^ = --J^Z = -S,ff (A.IO) 

P P deconfined 

Then using ( |A.9| ) and the approximations explained above the grand potential reduces to, 

72 °° 




P 

^ m=0 



(A.ll) 

+ ln(l 



In the limit that we are working in the sum in ( |A.ll ) receives primary contribution from 



large values of m, when the sum may be replaced by an integral. This approximation 
essentially corresponds to the fact that in the flat space approximation that we are working 
in the discrete spectrum on the sphere becomes continuous. Thus we have 



J dm m2|^ln(l -x"'+2+'^) +ln(l -x'^+s-'^n, (A.12) 
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Again neglecting the 3/2 in comparison to the large value of /i and performing an integration 
by parts we obtain 

which in the limit ^ 1 may be approximated by ^ 

n{f,,T) « J7o(m) - ^^rV', (A.15) 

Now to calculate the specific heat we use the following thermodynamic relations 

S = -{dTn)y^^, N = -{d^n)y^^, a = {dTE)y^^ = T{dTS)y^^ (A.16) 

Note that we hold the volume of the space {S'^) fixed everywhere and therefore it does not 
appear explicitly anywhere and we need not worry about it separately. Also it follows from 
the above relations that holding N fixed is the same as holding n fixed upto the leading 
order; therefore to leading order in fi/T we have 

a = 'Jf^T,^ (A.17) 

Note the linear dependence of specific heat on the temperature T which is a standard well 
known result. Also note that the specific heat is proportional to N"^. Since each fermi- 
surface has the chemical potential of order /x and contributes 0{T) to the specific heat, 
this implies that there are N'^ number of fermi-surfaces for the theory we consider here. 
A straightforward generalization of this argument to arbitrary dimensions shows that the 
main result of this appendix (presence of N'^ number of fermi surfaces) continues to hold 
in the generalized case. 

B. Scaling Argument 



We start with the Maxwell action ( p.lD coupled to the energy momentum tensor ( |2.6D of 



Fermi liquids with the profile ( ^.13 ) and ( 2.14 ). We define the metric with the AdS radius 



L and the gauge potential by the ansatz 

ds'^ = {-f{z)dt^ + giz)dz^ + k{z)de^ + z'^idx'^ + dy^)) , 

A = —h{z)dt. (B.l) 

We defined the rescaled quantities /3,m,p and p as follows 

p[z) = —^, p{z) = -^, l3 = —-f], m = -m. (B.2) 



^Here we have used the fact that in the /i/T 2> 1 hmit we have 
Jo e^r- -f 1 Jo D 



(A.14) 



-33- 



If we substitute (|B.2| ) into the equations of motion, the dependence on e, k and L drops out. 
Therefore in the bulk of this paper, we can set e = k = L = 1 without losing generality. 
If we would like to know the results for general values of e, k and L, then we just need 
to replace (f3,m,p, p) in the context of this paper with (/3,m,p, p). It is clear from this 
analysis that the normalized parameters m and (3 are not only proportional to the mass 
and species of bulk fermion but also depend on the charge e of the bulk fermion. 



C. The functions in the perturbative soliton star solution 

The functions constituting the linear order corrections in /3 to the solution inside the star 

f = ^{n^hl {hi + 8) -16) log {I -z') 

- 3hl {hl{8z + TT - 8 + log(64)) - Stt + 81og(64)) + 4 (3/i^ {hi - 8) - 16) tan-i(z) 
+ 4 {3hl {hi + 8) - 16) tanh"^(z) + 16(7r + log(64))^ 

9ln\^) = ^QQ^2^l,_^ f ( ~ 10/^" + ^) - + 1) + 2l0g(^ + 1)) 

+ (3/i§z(6(4z - 4 + log(32)) + Stt) - A0hl{z{~8 + Stt - 18 log(2)) + 8) 
-80/ioz(7r + log(64)) - 256z) + 20ho {-3h^ + 24/i^ + 16) z*tSin-\z) + 256^ 

= ( - (48z5 + (5^ - 52)z4 + 40z + 5n - 36) (C.l) 

+ 40hl ((4 + 37r)z4 - 8z^ + 37r + 4) 

+ lO/io {z" + 1) (^3 {hi + 8) /ig (^log Q (z2 + 1)^ + 21og(z + 1)^ 

+ {6hl {hi - 8) - 32) tan~i(z) + 8 (-2 log {z^ + l) - 41og(z + 1) + log(64))) 

+ 80/io ((4 + tt)z^ + 8z + tt - 12) - 512 (z^ - l) j 

= ^ ( {hi + 8) - 16) z (-4tanh-i (z^) + 4tanh-i(z) - log(4)) 

+ (tt (3/ig (/ig - 8) - 16) - 128) z + 4 {-3h^ + 2Ahl + 16) ztan^^z) + 128^ 

The 0{I3) correction to the star radius is given by, 
z(i) = --^ (-8(/io - m) {3hl - m^) + ^ {hi + l) (3/i^ - 6hlm^ - m*) 



-4 + 1) (3/^4 _ 6/1^^2 _ ^^^-1 l^py 
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The integration constants in the solution outside (^^) is given by 

= ^/io(^2(m-/io)'(3/io + m)log(^^-l^ -3/i^(-8 + ^ + log(64)) 

- 2Ahlm + Qhlm^{TT - 6 log(2)) + 2m^ (m^ - e/ig) (- log (/ig + m^) - 2 log(/io + m) 



+4 log(;io)) + log I + + ^ ) j + 4 {iht &hW - m^) tan-i 



/in 



+ m'^(7r + log(64)) 

1 3 
F2 = -Qho{ho - m) {iho + to) 



isTi ( ~ 3/ig(-39 + Stt + 30 log(2)) - 120/i^to + 10/ii^m2(l + Stt - 18 log(2)) 

+ 20 {ihl ~ Qhlrr? - hom^) tan"! + lOfto (^(3/i;] + Qhlw? - to^) (^log 



+2 log ( j j - (/lo - m)^(3/io + m) log (^-^ - 1 



+ 5/ioTO'^(-3 + TT + log(64)) + 8m^ 

i(ft.o - mf [mhl + 47/^o"^ + 24m2) 
" 1440 



Hi = -^(/lo - m)3(3/io + to) 

1 / / 2/1^ 

i/a = — 37r/i;^ - eTr/i^TO^ - 2(/io - to)^ (3/io + m) lo~ ' ° 



+ 4 (-3/i;5 + 6/i^to2 + TO*) tan^i (^^^ + {ih^ + 6/i^TO^ - to*) {-Ata.n\r^ 
+4tanh"* ( ) - log(4) ) + ShonT' - (8 + 7r)m* 



(C.3) 



D. Constants in the analytical solution near the tip of the soliton star 

In this appendix we list the constants in ( |5.1| ) as determined from the equations of motion 
The constant appearing in /(z) function 



_ 8 (-9/3/1^ + 20/3/i|]m2 - 15/3 /lom^ + 4 (/Sm^ - 45)) 
~ 9/3/i^ - lQ(5hlm? - 15/3/iom4 + 16 (/Jm^ - 45) ' 
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The constants appearing in g{z) are 

180 

^^"^^ ~ ~9(3hl, - 10/3/i|]m2 - 15(3hom^ + 16 (/Jm^ - 45) ' 

90 = — T ( - 1134/32C 

(9/3/i^ - 10/3/ii]m2 - 15/3/iom4 + 16 (/Jm^ - 45)) V 

+ 135/3^/1^ (26m2 + 15) - {866m^ + 1215) + 1188/3/i^ (/3m^ - 45) ^^''^^ 

+ {Um'^ + 81) - TeO/J/ij^m^ (/3m^ - 45) + 225/3^/igm^ (8m2 - 9) 

- 2220/3/10 m'' (/3m^ - 45) + 896 (/3m^ - 45)^^ 
The constant in k{z) is 

720 / 

k2 = o P (162/3C - 225;5/i^ (2m'^ + S) 

(/3 (9/i^ - 10/i|]m2 - 15/iom4 + 16m5) - 720) V 

+5/3/i[jm2 (I82m2 + 405) - 684/i^ (/3m^ - 45) - 75/3/i^m^ (lOm^ + 27) 

+1000/ii]m2 (/3m^ - 45) + 675/3/igm^ - 60/iom^ (/3m^ - 45) - 128m^ (/3m^ - 90)) 



259200 J 



(D.3) 



And the constant in h{z) is 

9/3/1^ - 10;5/igm2 - 15/3 /lom^ + 16 (/Jm^ - 45) 
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